Discovery of a Low-mass Companion to a Metal-rich F Star with the MARVELS Pilot Project

We report the discovery of a low-mass companion orbiting the metal-rich, main sequence F star TYC 2949-00557-1 during the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS) pilot project. The host star has an effective temperature T_(eff) = 6135 ± 40 K, logg = 4.4 ± 0.1, and [Fe/H] = 0.32 ± 0.01, indicating a mass of M_⊙ = 1.25 ± 0.09 M_⊙ and R = 1.15 ± 0.15 R_⊙. The companion has an orbital period of 5.69449 ± 0.00023 days and straddles the hydrogen burning limit with a minimum mass of 64 M_J , and thus may be an example of the rare class of brown dwarfs orbiting at distances comparable to those of "Hot Jupiters." We present relative photometry that demonstrates that the host star is photometrically stable at the few millimagnitude level on time scales of hours to years, and rules out transits for a companion of radius ≳ 0.8 R_J at the 95% confidence level. Tidal analysis of the system suggests that the star and companion are likely in a double synchronous state where both rotational and orbital synchronization have been achieved. This is the first low-mass companion detected with a multi-object, dispersed, fixed-delay interferometer

Theory of Dispersed Fixed-Delay Interferometry for Radial Velocity Exoplanet Searches

The dispersed fixed-delay interferometer (DFDI) represents a new instrument concept for high-precision radial velocity (RV) surveys for extrasolar planets. A combination of Michelson interferometer and medium-resolution spectrograph, it has the potential for performing multi-object surveys, where most previous RV techniques have been limited to observing only one target at a time. Because of the large sample of extrasolar planets needed to better understand planetary formation, evolution, and prevalence, this new technique represents a logical next step in instrumentation for RV extrasolar planet searches, and has been proven with the single-object Exoplanet Tracker (ET) at Kitt Peak National Observatory, and the multi-object W. M. Keck/MARVELS Exoplanet Tracker at Apache Point Observatory. The development of the ET instruments has necessitated fleshing out a detailed understanding of the physical principles of the DFDI technique. Here we summarize the fundamental theoretical material needed to understand the technique and provide an overview of the physics underlying the instrument's working. We also derive some useful analytical formulae that can be used to estimate the level of various sources of error generic to the technique, such as photon shot noise when using a fiducial reference spectrum, contamination by secondary spectra (e.g., crowded sources, spectroscopic binaries, or moonlight contamination), residual interferometer comb, and reference cross-talk error. Following this, we show that the use of a traditional gas absorption fiducial reference with a DFDI can incur significant systematic errors that must be taken into account at the precision levels required to detect extrasolar planets.Comment: 58 pages, 11 figures, 1 table, 3 appendices. Accepted for publication in ApJS. Minor typographical corrections; update to acknowledgment

Measurement of Spin-orbit Misalignment and Nodal Precession for the Planet around Pre-main-sequence Star PTFO 8-8695 from Gravity Darkening

PTFO 8-8695b represents the first transiting exoplanet candidate orbiting a pre-main-sequence star (van Eyken et al. 2012, ApJ, 755, 42). We find that the unusual lightcurve shapes of PTFO 8-8695 can be explained by transits of a planet across an oblate, gravity-darkened stellar disk. We develop a theoretical framework for understanding precession of a planetary orbit's ascending node for the case when the stellar rotational angular momentum and the planetary orbital angular momentum are comparable in magnitude. We then implement those ideas to simultaneously and self-consistently fit two separate lightcurves observed in 2009 December and 2010 December. Our two self-consistent fits yield Mp = 3.0 M_Jup and Mp = 3.6 M_Jup for assumed stellar masses of M* = 0.34 M_☉ and M* = 0.44 M_☉ respectively. The two fits have precession periods of 293 days and 581 days. These mass determinations (consistent with previous upper limits) along with the strength of the gravity-darkened precessing model together validate PTFO 8-8695b as just the second hot Jupiter known to orbit an M-dwarf. Our fits show a high degree of spin-orbit misalignment in the PTFO 8-8695 system: 69° ± 2° or 73°.1 ± 0°.5, in the two cases. The large misalignment is consistent with the hypothesis that planets become hot Jupiters with random orbital plane alignments early in a system's lifetime. We predict that as a result of the highly misaligned, precessing system, the transits should disappear for months at a time over the course of the system's precession period. The precessing, gravity-darkened model also predicts other observable effects: changing orbit inclination that could be detected by radial velocity observations, changing stellar inclination that would manifest as varying vsin i, changing projected spin-orbit alignment that could be seen by the Rossiter–McLaughlin effect, changing transit shapes over the course of the precession, and differing lightcurves as a function of wavelength. Our measured planet radii of 1.64 R_Jup and 1.68 R_Jup in each case are consistent with a young, hydrogen-dominated planet that results from a "hot-start" formation mechanism

Characterizing the Variability of Stars with Early-release Kepler Data

We present a variability analysis of the early-release first quarter of data publicly released by the Kepler project. Using the stellar parameters from the Kepler Input Catalog, we have separated the sample into 129,000 dwarfs and 17,000 giants and further sub-divided the luminosity classes into temperature bins corresponding approximately to the spectral classes A, F, G, K, and M. Utilizing the inherent sampling and time baseline of the public data set (30 minute sampling and 33.5 day baseline), we have explored the variability of the stellar sample. The overall variability rate of the dwarfs is 25% for the entire sample, but can reach 100% for the brightest groups of stars in the sample. G dwarfs are found to be the most stable with a dispersion floor of σ ~ 0.04 mmag. At the precision of Kepler, >95% of the giant stars are variable with a noise floor of ~0.1 mmag, 0.3 mmag, and 10 mmag for the G giants, K giants, and M giants, respectively. The photometric dispersion of the giants is consistent with acoustic variations of the photosphere; the photometrically derived predicted radial velocity distribution for the K giants is in agreement with the measured radial velocity distribution. We have also briefly explored the variability fraction as a function of data set baseline (1-33 days), at the native 30 minute sampling of the public Kepler data. To within the limitations of the data, we find that the overall variability fractions increase as the data set baseline is increased from 1 day to 33 days, in particular for the most variable stars. The lower mass M dwarf, K dwarf, and G dwarf stars increase their variability more significantly than the higher mass F dwarf and A dwarf stars as the time baseline is increased, indicating that the variability of the lower mass stars is mostly characterized by timescales of weeks while the variability of the higher mass stars is mostly characterized by timescales of days. A study of the distribution of the variability as a function of galactic latitude suggests that sources closer to the galactic plane are more variable. This may be the result of sampling differing populations (i.e., ages) as a function of latitude or may be the result of higher background contamination that is inflating the variability fractions at lower latitudes. A comparison of the M dwarf statistics to the variability of 29 known bright M dwarfs indicates that the M dwarfs are primarily variable on timescales of weeks or longer presumably dominated by spots and binarity. On shorter timescales of hours, which are relevant for planetary transit detection, the stars are significantly less variable, with ~80% having 12 hr dispersions of 0.5 mmag or less

Measuring Stellar Radial Velocities with a Dispersed Fixed-Delay Interferometer

We demonstrate the ability to measure precise stellar barycentric radial velocities with the dispersed fixed-delay interferometer technique using the Exoplanet Tracker (ET), an instrument primarily designed for precision differential Doppler velocity measurements using this technique. Our barycentric radial velocities, derived from observations taken at the KPNO 2.1 meter telescope, differ from those of Nidever et al. by 0.047 km/s (rms) when simultaneous iodine calibration is used, and by 0.120 km/s (rms) without simultaneous iodine calibration. Our results effectively show that a Michelson interferometer coupled to a spectrograph allows precise measurements of barycentric radial velocities even at a modest spectral resolution of R ~ 5100. A multi-object version of the ET instrument capable of observing ~500 stars per night is being used at the Sloan 2.5 m telescope at Apache Point Observatory for the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS), a wide-field radial velocity survey for extrasolar planets around TYCHO-2 stars in the magnitude range 7.6<V<12. In addition to precise differential velocities, this survey will also yield precise barycentric radial velocities for many thousands of stars using the data analysis techniques reported here. Such a large kinematic survey at high velocity precision will be useful in identifying the signature of accretion events in the Milky Way and understanding local stellar kinematics in addition to discovering exoplanets, brown dwarfs and spectroscopic binaries.Comment: 9 pages, 4 figures. Accepted for publication in Ap

Measuring Stellar Radial Velocities with a Dispersed Fixed-Delay Interferometer

We demonstrate the ability to measure precise stellar barycentric radial velocities with the dispersed fixed-delay interferometer technique using the Exoplanet Tracker (ET), an instrument primarily designed for precision differential Doppler velocity measurements using this technique. Our barycentric radial velocities, derived from observations taken at the KPNO 2.1 meter telescope, differ from those of Nidever et al. by 0.047 km/s (rms) when simultaneous iodine calibration is used, and by 0.120 km/s (rms) without simultaneous iodine calibration. Our results effectively show that a Michelson interferometer coupled to a spectrograph allows precise measurements of barycentric radial velocities even at a modest spectral resolution of R ~ 5100. A multi-object version of the ET instrument capable of observing ~500 stars per night is being used at the Sloan 2.5 m telescope at Apache Point Observatory for the Multi-object APO Radial Velocity Exoplanet Large-area Survey (MARVELS), a wide-field radial velocity survey for extrasolar planets around TYCHO-2 stars in the magnitude range 7.6<V<12. In addition to precise differential velocities, this survey will also yield precise barycentric radial velocities for many thousands of stars using the data analysis techniques reported here. Such a large kinematic survey at high velocity precision will be useful in identifying the signature of accretion events in the Milky Way and understanding local stellar kinematics in addition to discovering exoplanets, brown dwarfs and spectroscopic binaries.Comment: 9 pages, 4 figures. Accepted for publication in Ap

Application of the trend filtering algorithm for photometric time series data

Detecting transient light curves (e.g., transiting planets) requires high-precision data, and thus it is important to effectively filter systematic trends affecting ground-based wide-field surveys. We apply an implementation of the Trend Filtering Algorithm (TFA) to the 2MASS calibration catalog and select Palomar Transient Factory (PTF) photometric time series data. TFA is successful at reducing the overall dispersion of light curves, however, it may over-filter intrinsic variables and increase “instantaneous” dispersion when a template set is not judiciously chosen. In an attempt to rectify these issues we modify the original TFA from the literature by including measurement uncertainties in its computation, including ancillary data correlated with noise, and algorithmically selecting a template set using clustering algorithms as suggested by various authors. This approach may be particularly useful for appropriately accounting for variable photometric precision surveys and/or combined data sets. In summary, our contributions are to provide a MATLAB software implementation of TFA and a number of modifications tested on synthetics and real data, summarize the performance of TFA and various modifications on real groundbased data sets (2MASS and PTF), and assess the efficacy of TFA and modifications using synthetic light curve tests consisting of transiting and sinusoidal variables. While the transiting variables test indicates that these modifications confer no advantage to transit detection, the sinusoidal variables test indicates potential improvements in detection accuracy

Application of the Trend Filtering Algorithm for Photometric Time Series Data

Detecting transient light curves (e.g., transiting planets) requires high-precision data, and thus it is important to effectively filter systematic trends affecting ground-based wide-field surveys. We apply an implementation of the Trend Filtering Algorithm (TFA) to the 2MASS calibration catalog and select Palomar Transient Factory (PTF) photometric time series data. TFA is successful at reducing the overall dispersion of light curves, however, it may over-filter intrinsic variables and increase "instantaneous" dispersion when a template set is not judiciously chosen. In an attempt to rectify these issues we modify the original TFA from the literature by including measurement uncertainties in its computation, including ancillary data correlated with noise, and algorithmically selecting a template set using clustering algorithms as suggested by various authors. This approach may be particularly useful for appropriately accounting for variable photometric precision surveys and/or combined data sets. In summary, our contributions are to provide a MATLAB software implementation of TFA and a number of modifications tested on synthetics and real data, summarize the performance of TFA and various modifications on real ground-based data sets (2MASS and PTF), and assess the efficacy of TFA and modifications using synthetic light curve tests consisting of transiting and sinusoidal variables. While the transiting variables test indicates that these modifications confer no advantage to transit detection, the sinusoidal variables test indicates potential improvements in detection accuracy